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Abstract 

In this paper we present the results from a calculation of the full electroweak one- 
loop corrections for 7Z vector-boson pair production at hadron colliders. The cases 
of proton-antiproton as well as proton-proton collisions, at the Tevatron and the 
LHC, respectively, are considered. Results are presented for the distribution of the 
7Z invariant mass and for the transverse momentum of the final-state photon. The 
higher-order electroweak effects are numerically significant, in particular for probing 
possible anomalous gauge-boson couplings. 



1 Introduction 



Experiments at e + e~ and hadron colliders during the last decade have confirmed the pre- 
dictions of the Standard Model to very high accuracy (see for a recent review). In 
particular, the interaction of gauge bosons with fermions has been tested to a precision at 
the level of 0.1%. On the other hand, one of the most direct consequences of the Stan- 
dard Model as a non-Abelian gauge theory, the self couplings of the gauge bosons, are 
known with lower experimental precision. At the LHC, a substantial improvement in the 
measurement of the vector-boson self couplings will become feasible via the production 
of gauge-boson pairs providing the best direct tests of the non-Abelian gauge symmetry. 
Precise predictions from the Standard Model will allow to isolate possible deviations in 
the experimental data which may come from anomalous contributions to the gauge self- 
interactions indicating the presence of new physics beyond the Standard Model. If no 
deviations are observed, severe bounds on anomalous vector-boson couplings will be ob- 
tained (see e.g. [2] and references therein). Production of neutral vector-boson pairs, like 
7Z production, is well suited to search for non-zero ZZ'j and Z77 couplings, which vanish 
in the Standard Model, and could thus provide a clean signal of new physics. 

Previous studies of 7Z pair production in hadronic collisions were first performed at the 
level of the Born approximation based on the parton processes qq — > 7Z [3] and were later 
improved by the NLO QCD corrections [4j. The 0(a s ) QCD corrections to the hadronic 
cross section have been shown to modify the Born result considerably. These corrections 
are of special importance for high values of the 7Z invariant mass and for high transverse 
momenta of the photon, i.e. in the regions where also effects from possible anomalous 
vector-boson self couplings are expected to be identified as deviations from the Standard 
Model predictions. It is therefore important to have these predictions under control, which 
makes the inclusion of the higher-order contributions a necessity. 

Besides the QCD corrections, electroweak higher-order effects can become significant 
and have to be taken into account as well. The one-loop logarithmic electroweak effects in 
WZ and W'y production have been studied in [5]. Here we consider the hadronic production 
of 7Z pairs and present a calculation of the full electroweak one-loop corrections to the 
parton processes qq — > 7Z and to the hadronic cross section and distributions. Applying 
appropriate kinematical cuts, numerical results are given for the invariant-mass distribution 
of '-/Z and for the photon transverse momentum, in PP collisions at the LHC as well as in 
PP collisions at the Tevatron. 

2 Parton processes 

At the partonic level, pair-production of neutral vector bosons in lowest order can only 
proceed via quark-ant iquark annihilation, qq — > 7Z (q = u,d,s,c,b). NLO electroweak 
corrections consist of the one-loop virtual contributions and the contributions from real- 
photon bremsstrahlung . The subclass of loop diagrams involving virtual photons gives rise 
to infrared (IR) divergences; these are cancelled by including real-photon emission off the 
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Figure 1: One- loop diagrams with virtual photons (virtual QED contributions). The cross 
marks the QED mass-renormalization counter term from the virtual photon contribution 
to the quark self-energy. 

quark lines and integrating over the photon phase space. In the NLO-order electroweak 
contribution to the partonic cross section, da^ , it is convenient to separate the QED 
corrections from the residual part of the purely weak corrections, 

d4 } = + d 4 ED , (1) 

with the QED part split according to 

da^ ED = da^ irt + daj™ al . (2) 

Such a separation like in (0) appears natural for neutral-current processes since the virtual 
QED corrections form a gauge-invariant subset of the full one-loop diagrams; they are 
defined as the sum of the diagrams displayed in Fig. ^ Moreover, they are UV-finite 
when quark wave-function renormalization and mass renormalization, both exclusively 
from virtual photons, are taken into account, as indicated in the third line of Fig. 
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The QED corrections still contain mass-singular large logarithms of the type log s/m 2 q 
(with the invariant mass squared s of the parton process) arising from real and virtual 
photons collinear to a quark or antiquark with mass m q . These singular terms are universal 
and can be absorbed by a redefinition of the parton distributions functions (PDFs) of the 
initial-state quarks. By this procedure, the mass singularities disappear from the observable 
cross section, and the renormalised PDFs become dependent on the factorization scale 
controlled by the Gribov-Lipatov-Altarelli-Parisi (GLAP) equations [Hj. Those universal 
photonic corrections can be taken into account by a modifcation of the GLAP equations 
of QCD introducing an additional QED-evolution term [Zj. This leads to small corrections 
to the PDFs, at the per-mill level (2J 

The consistent treatment of the QED corrections would also require the inclusion of 
QED corrections in all the data used for global fitting of the PDFs. Current determina- 
tions of PDFs [HI E] do not include QED corrections, inducing thus an uncertainty which, 
however, should be small compared to the present uncertainties on the PDFs. 

Absorbing the collinear singularities associated with initial-state photon radiation into 
the PDFs introduces a QED factorization-scheme dependence. Our calculation is based on 
an explicit diagrammatic evaluation, with the collinear singularities factorized according to 
the MS scheme [TUj , but QED corrections are not taken into account in the PDF evolution. 

For the residual class of purely weak corrections, one only encounters loop diagrams 
and counter terms with massive virtual particles. Note that also in the counter terms 
for quark mass and wave-function renormalization all virtual photon contributions have 
been separated off as part of the QED corrections. The counter terms are determined by 
the electroweak renormalization scheme, which we choose as the on-shell scheme [TT| IT2] . 
The partonic cross section derived from the renormalized set of weak loop diagrams is 
hence UV-finite and does not contain any collinear or infrared divergences. The practical 
computation and evaluation of the complete one-loop terms was done with support of the 
packages FeynArts and FormCalc > where the version of the on-shell renormalization 
scheme is implemented. 

3 Hadronic observables 

The observable hadronic cross section for P + P(P) — > 7 + Z + X, with a given total 
hadronic CMS energy y/S, can be written as a convolution of the parton cross sections 
with the corresponding parton luminosities and summation over the various parton species. 
The one-loop order partonic cross section for qq annihilation is obtained from the lowest- 
order cross section dagJ and the NLO contribution (JTJ) after factorization of the collinear 
photon contributions in the QED part at the scale //p, 

da { ° +1 \s, t) = da^is, t) + da™* k (s, t) + daf D (s, t, p F ) , (3) 

with the invariant kinematical variables s for the partonic CMS energy squared and t for 
the momentum transfer between q and 7. Integration over t (applying appropriate cuts as 
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described below) yields the hadronic cross section in the following way 

dr 



1 r\r AB 
" Z (S)= I dr £ ^-a^\rS) (4) 

uu,dd,. 



°AB 



with the parton luminosity 

c\t r T T 1 

(5) 



q A (x,fj jF )q B {-,fiF) + qA(x,fJ jF )qB(-,^F) 
dr J T x 1 x x 

where qA(x,fip) [qa{ x i /4f)] denote the density functions for the quarks [anti-quarks] in the 
hadron A carrying a fraction x of the hadron momentum at the scale fip] {A, B) = (P, P) 
for the LHC and (P, P) for the Tevatron. The lower bound of the r-integration (r ) 
determines the minimal invariant mass of the parton system, s Q = tqS. In our case, To 
depends on the kinematical cuts applied. In order to have sufficiently large transverse 
poton momenta not too close to the beam axis, the following cuts 

Pi > pY, \y"\ < y 7,c > (6) 

are imposed for the photon transverse momentum pi^ and for the pseudo-rapidity y 1 of the 
photon, defined by 

Q 

V 1 = -log(tan-) , (7) 

where 9 is the scattering angle in the laboratory frame. The cut implies an energy and 
angular cut in the parton CMS, 

\ 2 

,T' C _L . /(^' c \2 I A/T2 



s> r S = W? + M 



cos,|<Wl-- 4 ^ C)2 



Mf) 2 • 

In the laboratory frame, |y 7 | can still become quite large. Then, the additional %f< cut 
restricts photons from coming too close to the beam axis in the laboratory frame. Assign- 
ment of specific values for the minimum j)\ and maximum \y c \, as used in our numerical 
analysis, are contained in Tab. Q 

Observables of particular interest are the distribution for the invariant mass M^ v = y/I 
of the j—Z final state configuration, 

dM7 z - s ^ ~dr~ °™ {TS) ' (9) 

u mv q t q 

and the distribution for the photon transverse momenta, 

r 1 dC AB Hrr {0+1) 

't jOs-af-V***- (10) 
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(TeV) 


pY (GeV) 




M H (GeV) 


LHC 


14 


100 


2.4 


115 


Tevatron 


1.8 


10 


2.4 


115 



Table 1: Input parameters for the kinematical variables and the Higgs boson mass entering 
the numerical evaluation. 



The kinematical constraints arising from the cuts divide the phase space into several 
subvolumes, over each of which is integrated separately. Integration over phase space 
was done numerically. To test the stability of the numerical integration, the results from 
an adaptive Gauss algorithm were checked versus those obtained from the Monte Carlo 
routine VEGAS7. As a further check, our results for the hadronic cross section in Born 
approximation were compared with the results of the earlier calculations in j^j. 

4 Numerical evaluation 

The results of the numerical analysis are of lowest order with respect to QCD in the 
parton processes. Hence, the numerical values are not yet the real predictions from the 
Standard Model; they are given here to point out and illustrate the effects of the higher- 
order electroweak contributions. The input parameters and the kinematical cuts have been 
chosen as listed in Tabled together with the factorization scale \Lp = 2Mz- 

QED and weak corrections were calculated separately. As already mentioned above 
for the QED corrections, the collinear singularities were factorized according to the MS 
scheme, but QED corrections are not taken into account in the quark distribution functions 
used for our analysis [8 . The residual non-singular QED corrections are quite small, at 
the level of 0.5%. There is, however, a left-over QED factorization-scale dependence, as a 
consequence of the missing QED evolution of the PDFs. This induces a scale uncertainty 
which is of the order of a few per mill as well. The calculation of the QED contributions 
thus provides more an estimate of the order of magnitude rather than precise numbers. 
From a practical point of view, the QED effects for 7Z production are negligible since they 
are small and covered up by the uncertainties of the quark distribution functions. 

The class of weak corrections induces more significant effects, in particular in those 
regions that are of interest for testing possible anomalous gauge couplings 0. Sensitive 
observables for the analysis of anomalous vector-boson couplings are the distributions of the 
hadronic cross sections in terms of the invariant mass of 7 and Z as well as the distribution 
of the transverse momentum of the final state photon. 

To illustrate the electroweak higher-order effects, we show in Fig. [21 the distribution of 
the 'jZ invariant mass in lowest-order approximation in comparison to the corresponding 
result including the full 1-loop electroweak corrections, both for the LHC and for the 
Tevatron. To be more quantitative, Fig. |3] diplays the purely weak corrections to the 
distribution, relative to the Born result. In the LHC case, the weak corrections reach a 
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considerable size around 20% for high invariant masses, whereas for the Tevatron they are 
in the lower percentage region. 

In a similar way, in Figures 0] and EJ the corresponding results are shown for the px 
distribution of the final-state photon. Again, the weak contributions are sizeable in case 
of the LHC, at the level of 30% and more for high p T . But also for the Tevatron case, they 
can exceed the level of 10%. The origin of the large effects at high invariant masses or 
high p T , respectively, are large logarithms, e.g. of the type a log 2 (s/Mf ), in the one-loop 
contributions to the parton processes. 

Finally, the electroweak loop terms depend also on the mass of the Higgs boson, which 
was kept at Mh = 115 GeV for the numerical results displayed in the figures. The de- 
pendence on Mh is not very strong, yielding a variation of a few per cent (relative) of the 
electroweak corrections setting Mh = 1 TeV. 

The electroweak corrections are negative whenever they are sizeable, opposite to the 
positive NLO QCD corrections [3], and thus partially compensate the NLO QCD contri- 
butions. 

To summarize, we have presented the results of a full electroweak one-loop calculation 
for the production of jZ pairs at hadron colliders, evaluated for both the Tevatron and the 
LHC. Like for other neutral-current processes, the electroweak 0(a) corrections naturally 
decompose into QED and weak contributions, which are separately gauge invariant and 
UV finite. The QED corrections, arising from virtual and real photon emission, contain 
collinear singularities which are factorized and absorbed in the PDFs. The remaining QED 
corrections are small, below the per-cent level, and insignificant in view of the uncertainties 
of the quark distribution functions. The weak corrections, however, are sizeable for large 
parton energies or transverse photon momenta, respectively, especially for the LHC, and 
have to be taken into account together with the NLO QCD corrections for testing the gauge 
structure of the Standard Model and probing the existence of anomalous gauge couplings. 



Acknowledgement 

W.H. thanks the Kavli Institut for Theoretical Physics, where this work was completed. 
This research was supported in part by the European Community's Human Potential Pro- 
gramme under contract HPRN-CT-2000-00149 "Physics at Colliders" and by the National 
Science Foundation under Grant No. PHY99-07949. 



6 



10" 



10" 



10" 



14 TeV 
100 GeV 
2.4 




> 

o 

O 



100 
50 

10 

5 



b 



f-.S 
T3 



Vs 



= 1.8 TeV 
= 10 GeV 
= 2.4 



500 1000 1500 2000 2500 100 200 300 400 500 

M7„v (GeV) Mf v (GeV) 

Figure 2: Distribution of the invariant mass in Born approximation (solid line) and 
including the full 1-loop electroweak corrections (dotted line), for the LHC (left case) and 
for the Tevatron (right case). 




Figure 3: Weak corrections to the distribution of the invariant mass relative to the 
Born result, for the LHC (left case) and for the Tevatron (right case). 
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Figure 4: Distribution of the transverse momentum in Born approximation (solid line) 
and including the full 1-loop electroweak corrections (dotted line) for the LHC (left case) 
and for the Tevatron (right case). 
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Figure 5: Weak corrections to the distribution of the transverse momentum relative to 
the Born result, for the LHC (left diagram) and for the Tevatron (right diagram). 
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